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Abstract 

We investigate quantum dynamical systems defined on a finite dimensional Hilbert space and subjected to 
an interaction with an environment. The rate of decoherence of initially pure states, measured by the increase 
of their von Neumann entropy, averaged over an ensemble of random pure states, is proved to be bounded from 
above by the partial entropy used to define the ALF dynamical entropy. The rate of decoherence induced by 
the sequence of the von Neumann projectors measurements is shown to be maximal, if the measurements are 
performed in a randomly chosen basis. The numerically observed linear increase of entropies is attributed to 
free-independence of the measured observable and the unitary dynamical map. 

PACS numbers: 03.65.Yz, 05.45.Mt 
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1 Introduction 

The notion of chaos in classical mechanics is well established, and any dynamical system characterized by positive 
Kolmogorov-Sinai entropy is called chaotic Q]- On the other hand it is not at all easy to generalize the definition 
of chaos for quantum theory E]. There exist numerous attempts to define the quantum counterpart of the 
Kolmogorov-Sinai entropy both for finite and infinite quantum systems (see @J and references therein). However, 
only two of them: CNT-entropy [5] and ALF-entropy |^j provide nonequivalent notions of quantum dynamical 
entropy which satisfy the following conditions: 

1) can be formulated in an abstract algebraic framework valid for general commutative (classical) and noncom- 
mutative (quantum) dynamical systems, 

2) coincide with the Kolmogorov-Sinai entropy when applied to classical systems, 

3) can be rigorously computed for several examples on noncommutative dynamical systems (different types of 
quantum shifts, automorphisms of quantum tori, fermionic quasi-free systems). 

In contrast to the coherent-states quantum entropy defined in 0] E] , both CNT and ALF entropies are always 
equal to zero for quantum systems with finite-dimensional Hilbert spaces. In particular this concerns the systems 
being quantizations of classically chaotic dynamical systems defined on a compact phase-space, we are going to 
analyze in this work. For the case of ALF entropy we can easily understand the mechanism leading to the apparent 
lack of "correspondence principle" for the K-S entropy. In fact we can see that classical-quantum correspondence 
exists provided the proper order of limit procedures is used. Recently this problem has been studied rigorously for 
both CNT and ALF entropies in the case of quantized Arnold cat maps [Hj- 

A typical quantum state coupled with an environment suffers decoherence, i.e. a generic pure state becomes 
mixed as a result of the non-unitary dynamics. As shown by Zurek and Paz [§], the initial rate of decoherence is 
governed by the classical dynamical entropy h. Vaguely speaking, any classical density evolving in a 2-dimensional 
phase space of a discrete invertible chaotic map T, is squeezed along the stable manifold and simultaneously stretched 
along the unstable manifold with the rate determined by the classical Lyapunov exponent A. In a similar way the 
corresponding quantum wave packet is stretched. So in a generic case, it becomes coupled with exponentially 
increasing number of states localized in the phase space. A natural assumption that these states are distinguishable 
by the environment implies the initially linear growth of the von Neumann entropy of a typical pure state, with the 
slope given by h[T]. 

Detailed investigation of the rate of decoherence in various setups is a subject of a considerable recent interest 
[Tfll fTTI Il2l 113] . The main aim of this work is to find a more precise bound for the increase of the average von 
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Neumann entropy in time. We analyze the decoherence in finite quantum chaotic systems subjected to the sequence 
of periodical measurement process and establish a link between the rate of decoherence and the partial entropy 
used for the definition of the ALF-entropy. 

The paper is organized as follows. In section [21 we recall the definition of the ALF quantum dynamical entropy. 
In section we discuss semiclassical limit of quantum maps and in section 0] analyze classical limit of the ALF 
entropy. In section we analyze the rate of decoherence and provide the upper bound (pTU) for the time evolution of 
the mean von Neumann entropy, averaged over the set of random initial pure states. These general results are used 
in sectionEUby studying the decoherence in a model system: periodically measured quantum baker map. Discussion 
of the results obtained in the context of the free-independent variables is provided in chapter 

2 ALF-dynamical entropy 

We describe an abstract, discrete-time reversible dynamical system in terms of the Hilbert space H, the unitary 
dynamical map U and the normalized vector |f2) € Ji. The relevant (generally complex) physical observables form 
a *-algebra A (i.e. a linear space with multiplication and adjoint operation) of bounded operators acting on H 
closed under the dynamics governed by U, i.e A 3 X UXU^ 1 = U(X) 6 A. We assume also that the state 
given by |0) restricted to the algebra A is time invariant, i.e. (Vl\U(X)\Vt) = (Cl\X\Q), In particular for the classical 
system (r, T, dj) with the reduced phase-space T, dynamical map 7 i— ► T(j), and the invariant probability measure 
dry, we obtain the Koopman's formalism with H — L 2 (Y,d^), (Uip)^) = ^(T^)), |0) — (function = 1) and A the 
algebra of all measurable bounded functions on the phase-space T treated as multiplication operators on L 2 (T,d , y). 
For a finite quantum mechanical principal system S the Hilbert space Ji = Hs <8> Ti-A describes the composite 
"system + ancilla" (S + A) with dimHs = dim^A = d, and the algebra A consists of elements of the product form 
X = Xs <8> 1a- The state \Q) is the "purification" of a generally mixed time invariant reference state of S. We have 
certain freedom in the choice of operator U. The simplest one is U = Us <8> 1, but it is sometimes convenient to put 
U = Us <8> Us, where Us is a properly defined "transposition" of Us such that U\fl) = \ This second choice is 
natural in the theory of infinite quantum dynamical systems usually formulated in terms of C* -algebras for which 
the present framework corresponds to the GNS-representation. 

ALF-entropy is defined in several steps using the notion of partition of unity X = {X\, X 2 , ■ ■ ■ , X^; Xj € A, 
XjXj = 1}. The partition of Hs generates a corresponding tensor-product partition of the composed space 
H, namely X (£> 1 = {Xi ® 1, . . . , Xk <8> 1}. To shorten the notation we will write in the sequel X instead of X ® 1. 
Partitions of unity can be composed, X o Y = {XjY m } and evolved U(X) — {UXjU^} to produce finer partitions, 
X 4 = W t_1 (X) o • • • oW(X) o X. We use the following notation for the multi-time correlation matrices given by 

^1,, Jt - {n\xlu{x} 2 ) . . .u^ix^u 1 - 1 ^) . . .u(x i2 )x h Q) 

= (n\xltfx} 2 tf . . . x]x lt . . . ux l2 ux n n) . (1) 

Here erpC'] is a positively defined, k l x fc* complex- valued matrix with a trace equals one. By 5t[X, U] we denote 
its von Neumann entropy 

S t [X,U] = -tr(o-[X*] IncrpC*]) = -tr^X^lnf^X*]), (2) 

where f2[X*] is a finite range density operator acting on H 

k 

n[X t }= J2 I UX i> ■ ■ ■ UX h ux n V) (UX jt . . . UX h UX n 0| . (3) 

h,h,— jt=l 

Equality follows from the fact that the spectrum of the operator J2j=i li)O I i s identical (including degeneracies) 
to the spectrum of the fc x k matrix [{i\j}] except the eigenvalues equal to zero. 

For any partition of unity X, one may introduce the corresponding dynamical map <I>x in the Schrodinger picture 

k 

P -> $*{p) =Y, X iP X l- (4) 

3=1 

This map sends an arbitrary density operator p in the set of density operators. Iterating the state |^)(^| i-times 
by the map ^uk. we obtain 

ft[x<] = [Wflnxni) , (5) 
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where UX = {UX X , UX 2 , UX k }. 

The formula of above suggests a new interpretation of St[X, U] as the entropy of the density matrix obtained 
by repeated measurements performed on the evolving system plus ancilla with the initial pure entangled state |f2). 

The dynamical entropy of partition h[X, U] is defined as a limit 

h[X, U] = limsup -S t [X, U] . (6) 

t— >oo t 

Finally, the dynamical entropy of U is a supremum over a given class {X} of physically admissible partitions [gj 

h[[/]=suph[X,[/] , (7) 

{X} 

i.e. over all generalized measurement processes. One can choose as {X} all partitions in the algebra A or restrict 
the supremum to bistochastic ones (satisfying ^XjX^ = 1), or even partitions with Xj proportional to unitary 
Uj, namely Xj — PjUj, Pj > 0, J2jPj = 1- Although in general one obtains different types of ALF-entropies they 
all coincide in the classical case and for all known examples of noncommutative dynamical systems. 
Since for any finite quantum system the entropy St[X, U] is limited 

S t [X,U] < min{ilnfc,d} , (8) 

the asymptotic rate © of the entropy production gives zero, independently of the investigated unitary dynamics 
U. On the other hand, one may analyze the initial rate of the entropy St, which at small times was shown J3j to 
be determined by the classical entropy h. 



3 Classical limit of quantum maps 

Consider a classical dynamical system with a compact (reduced) phase-space T equipped with the probability 
measure dy. This measure is assumed to be invariant with respect to the dynamical map 7 1— > T(/y) . We say that 
this system (T, T, dry) is a classical limit of the sequence of quantum systems if: 

a) there exists a sequence (C d( - n \U n ;n = 1,2,...) of d(n)~ dimensional Hilbert spaces and d x d- unitary 
matrices, 

b) there exist a quantization procedure which to any real function 7(7) (usually satisfying some smoothness 
properties) associates the sequence of self adjoint operators F^ 71 ' acting on C d . 

c) for any set of observables fx, f%, ■ ■ ■ , fk and any sequence of time steps ti, t%, . . . , tk the correlation functions 
converge, 

lim -^Tt(L#^ n >[r-«i[#^ = / d 7 f 1 (T t H7))f 2 (T t H7))---fk(T t Hl)) ■ (9) 

Both, the quantization procedure fj 1— > Fj and the choice of U n are not unique. The maximally mixed states 
= t/d correspond to the uniform normalized measure dy. 

The example of such structure has been rigorously studied for the Arnold cat maps in the recent paper :8|. 

To present an example of a family of quantum maps, corresponding to certain classical system, we are going 
to recall the construction of quantum baker map, originally due to Balazs and Voros 23 and modified later in 
^SlElEj- Classical baker map is defined as a transformation of a unit square - the compact phase space T with 
the coordinates q (position) and p (momentum), 

T3j=(q,p)^ T B (i) = (2q - [2 q ], (p+ [2g])/2) G V, (10) 

where [2q] denotes the integer part of 2q. This map is hyperbolic and its Kolmogorov-Sinai entropy is equal to 
In 2. Such a transformation may be quantized in a finite Hilbert space C d . In an ordered orthonormal basis named 
position eigenbasis {ej} we introduce a periodic translation operator 

Uej=e j+ i, j = 1, . . .,d- 1 , Ue d = e x . (11) 

Diagonalization of U leads to the conjugated basis - momentum eigenbasis {e^}, 

Ue k = exp(2mk/d)e k ■ (12) 
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Analogously to 11 111 the translation operator in the momentum eigenbasis is introduced, 

Vet = e k +i , Ved = ei , k = 1, . . . , d - 1 . (13) 
This operator is diagonal in the position eigenbasis, 

Vej = cxp(— 2mj/d)e ;, (14) 
and the transformation between position and momentum basis is given by the discrete Fourier transform 

h = E = E 7^ e ~ 2 ^' /d e ^ ■ ( 15 ) 

Having defined the group of translation operators corresponding to a classical torus it is possible to link 

the unitary operator 

&-(W-(*jjf* (16) 

acting on C d , where ci is even integer (e.g. d — 2n), to the classical transformation defined by iflOjl . The translation 
operators allow one also to define a finite dimensional operator corresponding to any classical observable described 
by a continuous function / on T. Let us defined the Fourier expansion of /, 

/( g ,p) = E a jk e~ 2 ^ q/d ^vld (17) 

Thus the operator corresponding to observable / may read as follows 

FW = E <hk V 3 U k . (18) 

j,k 

As it was mentioned above the quantization procedure is not unique. Another set of operators F^ n > may be obtained 
if we use different ordering of translation operators, since they do not commute, UV = VXJe llxl l d . It is also possible 
to generalize whole quantization procedure by introducing translation operators which are not exactly periodic, but 
periodic up to a phase factor (e.g. U d = e 27Tlx ^ d 1, V d = e 2mx ^ d 1) [TT)|IT8|. 

Different properties of such a quantum baker map were studied in ^3 QUI El EH , and the correspondence with 
the classical system (fflTTfl was established. Although we cannot provide a formal proof that the quantization iflfill 
satisfies the property (jHJl, we are going to use this model in further numerical investigations. 



4 Correspondence principle for dynamical entropy 

For any finite dimensional quantum system inequality (HJ holds, so the quantum dynamical entropy h[U) — 0. This 
fact is sometimes interpreted as the lack of correspondence principle for dynamical entropy. In section we defined 
a family of quantum maps, parametrized by an integer index n, such that in the semiclassical limit n i— » oo the 
dimension d(n) becomes infinite. Taking a sequence of quantum systems (C d , U n ) with the classical limit (r, T, dj) 
in the sense defined above ®, we may start with a functional partition of unity f = {ft, f%, ■ ■ ■ , fk) X^' \fj(l)\ 2 = 1} 

and construct its quantum counterparts F„ = {F^ 1 , Fj, , . . . , F^} using a suitable quantization procedure. The 
entropy St[F n , U n ] is computed using Eq. Q, 

tt[K] = [QUnlfJ ® l(|*n)<*n|) , (19) 

where 

1 d 

I**) = -f= d E l e ™) ® IO > {I4>} " baSiS 111 Cd ( 2 0) 

?n — 1 

is the purification of the tracial state 1/d of the system given in terms of the maximally entangled vector in C d ®C d . 
Then according to the Eq. © 

S t [f,T] = hm S t [f n ,U n ] , (21) 
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the classical dynamical entropy of the partition can be recovered by taking first the classical limit n — ► oo and then 
long time limit t — > oo. Usually, for a given classical system with the K-S entropy h[T] there exist many "optimal 
partitions" f with lnfc > h[T] for which St[t,T] « t ■ h[T] with a given accuracy or even exactly ( generating 
partitions, Markovian partitions |2|). Therefore we can expect that for large enough n and the optimal choice of the 
partition f the entropy S t [F n , U n ] displays linear growth with the rate given by the K-S entropy h[T] for t below 
tmax = 21nd/h[T] and then saturates at the maximal value 2 had. For the regular dynamics T with h[T] = we 
expect a slower (logarithmic) increase of St[F n , U n ] up to the maximal value. 

In formula I|19J1 any purification of the tracial state can be used while for concrete computations some choices 
could be better than the others. The set £ of pure states {\a) : (a\a) — 1} in C d equipped with the probability 
measure da and satisfying 

da\a)(a\ = -1 , (22) 
e « 

will be called complete set of vectors. The natural examples of £ are orthonormal basis or coherent states generated 
by the irreducible representations of certain compact Lie groups on C d . The distinguished example is the set of 
all pure states Vd = (DP d_1 with the natural unitary invariant probability measure (Fubini-Study measure). The 
following theorem provides the most general representation of maximally entangled vectors in terms of complete 
sets. 

Theorem 1 

1) Any purification of the tracial state t/d on C d is a maximally entangled state given by 

1 d 

l*> = -7= J2 l e ™> ® l e "> ' {l e ™»' - basis in Cd • (23) 

m— 1 

2) There is one to one correspondence between antiunitary matrices acting on C d and maximally entangled 
states in C d ® C d which can be expressed in terms of an arbitrary complete set of vectors £ as 

G <-> |* G ) = VdJ da\a) <g> \Ga) . (24) 

Proof 

1) Follows directly from the Schmidt decomposition. 

2) Take two complete sets £ , £' and two antiunitary matrices G,G'. Define \^g) by H24I1 put |^g') = 
Vd f £/ dj3\(3) <8) \G'0) and compute the scalar product 

(^ G ,\^ G )=d [ da j dp(0\a)(G(G- 1 G')P\Ga)= [ dp{(3\{G- 1 G')f3) . (25) 



From J2U it follows that ||*G'|I = ll*c?|| = 1 and *g' = *G if and only if G — G' . Obviously, for any vector of the 
form l(2"3jl we can choose the unique antiunitary matrix satisfying Gej = e' . 



5 Entropy production as a measure of decoherence 

For a generic quantum system S interacting with an environment (e.g. measuring apparatus) its initial pure state 
becomes mixed due to the increasing system-environment entanglement. Assuming that the reduced dynamics is 
given by a completely positive map Q, to measure the decoherence we may use the von Neumann entropy 

E[X,a] = 5(*x(|a)(a|)) = S(a a [X]) , (26) 

where |a) € Hs is an initial pure state of the system and <r Q [X] is k x k correlation matrix with element 
(a|XjXi|a). We extend this construction to the case of discrete time finite quantum dynamical system with the 
unitary evolution U interrupted by a measuring process (or generally interaction with an environment) described 
by the partition of unity X or equivalently by the 

E t {X,U,a] = S({<S> U x} t (\a)(a\)) = S( CT Q [X']) . (27) 

The quantity of above, bounded by ln(dim7is) can be strongly dependent on the initial state of the system. 
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Assume now that the system S is finite i.e. Hs = C . In order to obtain a more universal measure we can 
average the entropy over a complete set £ of pure states {\a}} . The entropy averaged with respect to £ is equal to 

Et [X, U,£} = J daE t [X, U, a] < In d , (28) 

and its increase (entropy production) characterizes the magnitude of the decoherence process. 

Since the entropy E t is bounded from above, its asymptotic production rate, E t /t tends to zero for t — > oo. On 
the other hand, we will be interested in the initial production rate. Studying a discrete dynamics we cannot define 
the derivative dEt/dt, but we may for instance study the entropy produced after each initial time step. Analyzing 
the trivial dynamics, U = 1, and measurement process governed by projection operators, Xj = Pj = {Pj) 2 , the 
entropy is produced only once, and E t = E\ for all t > 0. Therefore, to characterize in this situation the unitary 
dynamics, and not the measurement process itself, we are going to use the quantity AE = E2 — E\ as a measure of 
the decoherence. For comparison we define the initial production of the partial ALF-entropy, AS = S2 — S\. 

Defining the ALF entropy, which characterizes the unitary evolution U one uses the supremum over all 
operational partitions of unity. Let us emphasize that there is no point in performing such a step by studying the 
initial decoherence rate AE. Since the set of transformed operators, Pj — > P- — PjV (P = {Pi, . . . ,Pk}) with 
arbitrary unitary V is also a valid identity resolution, then E t [P,U,£] — E t [P' , VU, £] so the supremum over all 
possible measurements will be independent of the unitary dynamics U studied. 

It follows from Eq. (|2j and Eq. © that the time dependent entropy S*t[X, U], which appears in the definition 
of the ALF-entropy and in the semiclassical regime is related to the KS entropy, describes also the magnitude of 
a certain decoherence process. However, this process involves maximally entangled state of the system plus ancilla 
while the natural decoherence measure should be defined in terms of the system alone, like E t [X, U, a] or E t [X, U, £} . 

In order to compare both entropies E t [X-, U,£] and St[X, U] we need the following technical result. 

Take the dynamical map $y defined by the /c-elements partition of unity Y = {Yi, Y2, . ■ ■ , Y/.} as in Q and 
an arbitrary complete set of vectors £. We use the notation a a [Y] and <r[Y] for the k x k correlation (density) 
matrices with matrix elements (a|Y^y;|a) and ^Tr(Y^li) respectively. We introduce also the tracial norm \\A\\x = 

Tr[(^At) 1 /2] an d the Hilbert-Schmidt norm ||vl||2 = [Tr(AA^)] ^ 2 for a matrix or an operator A and the entropy 
function rj(x) — —a; In a;. 
Theorem 2 

A > S{a[Y]) - J daS(a a [Y]) > B , (29) 

where 

A= I d^HYl-^MIUlnfc + ?7(HY] -^[YJIK)) , (30) 

B = imax{^rf a (||a[Y]-a«[Y]|| 1 ) 2 ,^da(||a[Y]- ( 7 Q [Y]|| 2 ) 2 } . (31) 
Proof. We use the inequalities for the relative entropy |2"l l21l l3"] 

S(p\u)=Tv(p\np-p\nw) > ~max{||p - \\p - cj\\ 2 2 } . (32) 

Putting p = <J a \Y] and u> — a[Y] = J £ daa a [Y] and averaging over da we obtain the lower bound IpU- The upper 
bound ll-'iOII follows directly from the Fannes inequality pOd], 

\S(p)-S(u;)\<(\\p-u\\i]n{^n)+vQ\p-u;\\i) ■ (33) 

A basic consequence of the Theorem 2 is the inequality 

21nd> St[X® l d ,U] > E t [X.,U,£] < Ind , (34) 

and the fact that StpC, U] = E t [X-, U,£] if and only if <r Q [X*] = cr[X*] for almost all a (except of a set of measure 
zero). Numerical results show that for small times both quantities are comparable St[X, U] ~ £t[X, U,£]. We can 
provide some arguments in favour of this behavior in the case of the complete set Va of all pure states. 
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Consider the fluctuations of the matrix elements of the fc x fc matrix a a [Y] treated as random variables with 
respect of the uniform measure over the set of all pure states Vd- Deviation of a matrix element from its average is 
given by 

5=m* a [Y]-a[Y]\0\ . (35) 
where |£) is an arbitrary normalized vector from C fc . Expectation value of the operator cr Q [Y] is equal 



(£K[Y]|0 = <a|££&^l 



(36) 



The positive operator B = ^2ii^jY^Yj < 1, i.e. for any normalized vector \<f>) £ C d , {4>\B<p) < 1. So, it can be 
written in the form of convex sum of projectors into its eigenvectors i.e. B — J2i bi\^i) (0 < bi < 1). Let ai 
denote coefficients of the random normalized state \a) with respect to the eigenvectors of the operator B. Finally 



(37) 



The numbers (|a;| 2 ~ d^ 1 ) take positive and negative values of the order c? _1 but sum up to zero. However, when 
multiplied by another random variables bi € [0, 1] they behave like "steps of the random walk" yielding a sum of the 
order \fd x dr 1 and hence 

(38) 

Therefore the fluctuations of the norm ||<r[Y] — <j a [Y]||i behaves like k/yd . 

In the time-dependent case it means that for fc* -C d we have S*[X, U] ~ £?t[X, U, £]. Moreover, a random choice 
of |a) € V d gives typically S{a a [l^]) ~ E t [X, U,£] . 



6 Decoherence in periodically measured baker map 

To illustrate the results presented in previous section on a concrete example we analyze the quantum baker map iflfill 
subjected to periodic sequence of measurement performed in the momentum basis. Entire, non-unitary dynamics 
of the system is described by the superoperator 

P' = J2 P f U P ufp f ■ (39) 

i=i 

The set P p of fc projection operators fulfills the identity resolution, Y^j Pf — 1) since the measurement process 
corresponds to the partition of phase space into k equal intervals in momentum, 

{jd/k ~\ 
pf-pf= E 1^1 . ( 4 °) 
t=(i-l)d!/k+l J 

where ej are the momentum eigenstates defined by ltl2ll in C d , and the size d of the Hilbert space is an integer 
multiple of fc. 

Iterating numerically quantum map H39|l we compute how both entropies St[X, U] (SJ and E t [X., U,£] l|28|l vary 
in time. Fig. presents the initial growth of the entropy ^[P 1 ",?/] and E t [P p , U, Vd] , where we have averaged 
the entropy over the entire set of pure states with the natural measure, £ — Vd- As the evolution operator U we 
took the quantum baker map Ub defined by l(T6|l part (a), (c) and (e), and its square Ub 2 (b), (d) and (f). The 
partition P p is composed of fc = 2, 4, 8 projection operators. Entropy E t is averaged over a sample of 32 randomly 
chosen pure initial states. To guide the eye we plotted solid lines corresponding to growth with the rate of classical 
KS entropy, which is equal In 2 in the case (a), (c) and (e) and 2 In 2 in (b), (d) and (f). The slope of the dashed 
lines is equal to the maximal allowed growth of entropies, equal to lnfc. Aiming for the semicalssical regime, we 
have taken the maximal dimensionality of the Hilbert space, which was allowed by the computer resources at our 
disposal. In order to compute entropy S t one has to diagonalize matrices of size d 2 , so we could work with matrices 
size d = 64 = 2 6 . To obtain the entropy E t one needs to study the time evolution of density operators acting on 
ttd, so we succeed to work with systems of the size d = 512 = 2 9 . Entropy E t obtained for d = 64 is smaller than 
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St according to the analytical bound Il-Ull . As discussed in the size d of the Hilbert space determines only the 
saturation level (E t (t — > oo)), but does not influences the initial entropy rate. If the measurement scheme is tuned 
to the classical dynamics, i.e. case (a) and (d), the rate of the initial growth of both entropies coincides with the 
classical dynamical entropy h[T] of the map which is equal In 2 for the baker map, and 2 In 2 for its square. If the 
resolution of the measurement is not sufficient — (b) the classical chaos cannot fully manifest itself and AE and 
AS are smaller than h[T] and equal to Ink. In the opposite case, panel (c), (e) and (f), a finer resolution of the 
measurement (lnfc > h[T]) allows for the decoherence with the rate faster than it can be expected from the classical 
entropy. Hence such a measurement can be responsible for entropy production faster then it may be predicted 
basing on the degree of the classical chaos. As visible in panel (e) this effect is larger for the entropy S t . 




Figure 1: Initial growth of entropies St[P p ,/7] (o) and E t [P p , U, Vd] (x) computed for baker map U = Ub (a), c) 
and e)) and its square U = Ub 2 (b), d) and f)) where partition P p corresponds to division of classical phase space 
into k = 2(a) and b)), k = 4((c) and d)), k = 8((e) and f)) equal intervals in momentum coordinate. 

To demonstrate other features of the measurement process we investigated the time dependence of both entropies 
St[P R 7 U] and E t [P R , U, Vd] for different choices of the partitions P R . Fig. El shows the initial growth of both 
entropies calculated with the same evolution operators as in Fig. However, the partition P R was obtained by 
rotating the projective partition P p by a random unitary matrix V, namely = VPFV\ for all j = 1, . . . , k. 
The label R decorating the symbol P R of the partition emphasizes the fact that the measurement is performed in a 
random basis. Such a measurement will give k different results with equal probabilities tr(P^)/c? = 1/fc. In Fig- EJ 
both entropies initially increase with a nearly maximal slope which is equal to lnfc = In 8 (dashed line). Unitary 
evolution operator U does not influence the behavior of both entropies, and the data presented in both plots (a) 
and (b) hardly differ. 

The fact that the rotated partitions leads to an almost maximal allowed growth of both entropies may be 
explained by the following argument. Both St[X, U] and E t [X, U, Vd] may be expressed as the von Neumann 
entropy of a mixed state obtained by the operator [$[/x]* applied to a pure state (see Eqs. © and ffijl ). Let 
us denote the rotated partition by Y = V~KV\ where V is an arbitrary unitary matrix. From definition of the 
operator $x GJ we have 

[*inr]*(p) = [*irvxvt]*(p) = V&viuTrdKVi pV)V* = U[<Vx] (41) 

where the initial state p becomes rotated, p' = V^pV, and the evolution U is replaced by U' = VWV. The von 
Neumann entropy depends only on the spectrum of density matrix, so S([^uy] t (p)) — 5([$[/'x]'(p')- I n the case 
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Figure 2: As in the Fig. Q](e) and f) with k = 8 for the measurement in a random basis P R . 

of the entropy St the evolution operator is a tensor product U = Us ®Ua, so are the partition X = Xg £§) 1, and 
the matrix V = Vs <8> 1. Since we choose randomly only Vg, the state of the ancilla is not important. When p 
is a maximally entangled state so is p' , since those state are invariant under local operations. Hence we obtain 
S t [Y, U] = St[X, U']. The case of the entropy E t is simple, since here the averaging over all initial pure states \a) 
automatically cancels the difference between p and p' . Therefore the equality Et [Y, U, VJ\ = E t [X,U',Vd] holds. 
As we can see in both cases the randomly rotated partition Y is equivalent to model with original partition X 
and the evolution operator rotated into U' = V*UV, where V is a random unitary matrix, generated according to 
the Haar measure on U(N). Although the spectra of U and V are equal, the eigenvectors of U' are random, and 
the operator $(j'x generates the maximal growth of the von Neumann entropy AE w minjlnd, lnfc}. Here hid 
represents the maximal von Neumann entropy of the mixed state, while the measurement X transforms any pure 
state into a mixture with entropy bounded by In k. 

This argument shows that for finite systems taking the supremum over all possible partitions of unity leads to 
the maximal allowed growth of both entropies irrespectively of the analyzed unitary dynamics U. Note that the 
upper bound In d is only slightly larger that the average quantum dynamical entropy of a random unitary matrix 
U distributed according to the Haar measure on U(d) 0. 

On the other hand, one may pose a question, how to restrict the set of possible partitions, such that the rate 
of growth of von Neumann entropy could correspond to the dynamical entropy of the classical system. To analyze 
this problem compare the properties of the momentum partition P p and the random partition P R in the phase 
space. To represent the partition member we make use of the Husimi-like representation, 

x ] (q,p)=(q,p\X}X 3 \q,p). (42) 

Here \q,p) denotes the Gaussian states localized on torus, the same as in Ref. ^lEOI- In Fig. we show the 
phase-space representation of two partitions of unity P p = {Pi, P 2 }, and P R = {VP\V\ VP 2 V^}, each consisting 
oik — 2 operators. P p denotes the partition into equal intervals in momentum coordinates and P R is a partition 
into two subspaces of equal size determined by a random unitary matrix V, As may be seen in the picture, the 
operators Pf and P 2 P are by construction localized in lower (upper) region of the phase space, while Pf 1 and P 2 R 
are totally delocalized. 

These results suggest that in order to predict the decoherence rate AE (and AS) one should consider only these 
partitions P c , which have a well-defined classical limit. This is the case if each of the operator Xj is localized on 
a subset tj C T of phase space, so in the classical limit Xj(q,p) — > XeAQiP), where Xe,- denotes the characteristic 
function of €j . 

The discussion of above has some important consequences for the decoherence processes in quantum systems 
having chaotic classical limits with the K-S entropy h[T] > 0. Namely, if only the interaction with the environment 
can be described by a partition of unity having a well defined classical limit (and In A: > h[T]) we expect that 
the decoherence effects give the entropy production per single time step of the order h[T) for the generic initial 
conditions. The linear entropy increase has to break down for times t of the order of In d/h[T] On the contrary, 
if the partition of unity is chosen randomly and has no well-defined classical limit we do not expect any restrictions 
on entropy production rate, except the general upper bound AS < In k, related to the number k of the Kraus 
operators. 
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Figure 3: Phase-space representations of partitions P p = {Pf , -Rf} an d P R = {Pi > P^} consisting of k = 2 opera- 
tors. The real function Xj(q,p) — (q,p\XjXj\q,p) is plotted in dark scale, where X, denotes one of Pf , P£, P/ 1 , Pf. 
The partition P p corresponds to partition on upper and lower half in momentum coordinate while P R is partition 
corresponding to the projections in a a randomly selected basis. Note that the coherent states representations of 
P/* and P2 are delocalized in the phase space. 

7 Occurrence of free-independent variables 

Semiclassical arguments of the previous sections do not explain the striking phenomenon observed in numerical 
computations of St[X, £7] and Et [X, U, Vd] as presented in Fig. G] (a) and (d) and Fig. Namely for the two 
situations: 

a) the quantum system with chaotic classical limit and the semiclassical, projection valued choice of the partition 
P = {Pi, P 2 , • • .,Pfc},trPj = d/k satisfying Ink < h[T], 

b) nontrivial U and the random choice of the partition P. 

St[X, U] grows almost exactly linearly like tin A: and then rapidly saturates at the maximal value 2 hid. The 
entropy Ef [X, U, a] with a random choice of |a) follows the same plot up to its maximal value In d. This means 
that in both cases a) and b) the correlation density matrices possess a very special structure corresponding to the 
maximal admissible entropy, 

v[V% x ,...,ir,h,...,h = 'I '"•(/' ' T /'.,'' • ■ • Ph p h ■ ■ ■ ' /',' /' ) - ; . • • • S Ujt , (43) 
and similarly for a typical vector \a), 

a a \P% u ...^ u ..., Jt = (a\(P jl tfP h UK..P jt P it ...UP i2 UP il )\a) ~ ^5 lin ■ ■ ■ 5 Hn , (44) 

under the condition of fe* -C d. 

The simple form of the correlation functions l|43|) and 14411 suggests the existence of a certain statistical law 
satisfied by the noncommutative variables {Pj, U} with respect to the tracial state or a typical pure state |a). This 
law should be strictly obeyed in the limit d — ► 00 but even for relatively low dimensions reproduces the data with 
very good accuracy. Such situations are common in Nature. Gaussian and Poisson probability distribution are 
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very successful in describing experimental data while their rigorous derivations involve limit theorems with strong 
statistical independence assumptions. 

We advance the following statistical hypothesis: 

For both cases a) and b) and large Hilbert space dimensions d the operators {A, U} behave asymptotically like 
free-independent random variables with respect to the tracial state or a typical pure state. 
Here A is an arbitrary observable with the spectral measure (Pi, P2, . . -Pk)- 

We have to explain now the notion of free-independence. In the classical probability theory (complex) random 
variables form a commutative *-algebra and the probability measure defines a positive normalized functional / 1— > (/) 
- the average value. The random variables f x, f 21 ■■■ 1 fn are called statistically independent if 

(/i/ 2 .../„) = (/i)(/ 2 ) •■•</„) • (45) 

In noncommutative probability the basic object is a unital generally noncommutative *-algebra A with a state 
(positive normalized functional) 4>. Due to the noncommutativity the average 4>(xiX2 ■ ■ ■i m );ij € A depends on 
the order of random variables. Hence, the direct extension of the definition 14511 is not very interesting and essentially 
corresponds to product states on tensor product algebras. Instead, in noncommutative probability we have different 
notions of independence which take into account possible algebraic relations between random variables (e.g. CCR, 
CAR, etc (221 )■ I n the last decade the so-called free families of random variables (or free-independence) introduced 
by Voiculescu (23 E3| attracted attention of physicists mainly due to the relations with random matrices theory. 

Denote by w(x) an arbitrary polynomial in variables x, x^ E A. The collection of noncommutative random 
variables x\, x%, ■ ■ ■ , Xu is called free-independent if 

<j){wi{x p ( 1) )w2{x p ( 2 )) ■ ■ ■ w m {x p(m) ) =0 , (46) 

whenever t/>(uij (x p{]) ) = and p(j) ^ p(j + 1) for all j = 1, 2, . . . , m, p(j) G {1,2,..., A:}. 

It has been proved that the Wigner semicircular probability distribution is a consequence of the central limit 
theorem for free-independent random variables similarly to the origin of Gaussian probability distribution in the 
context of statistically independent commutative variables. Moreover, the free-independent variables naturally arise 
as limits of large random matrices |23l 124] . The consequences of free independence are illustrated by the following 
example. 

Example. Take a family of orthogonal projections Pi, P2, ■ ■ ■ , Pk and a unitary U, all from the algebra A with 
the state 4>. Assume that for any A = ^2 a jPj the pair of random variables {A, U} is free-independent and moreover 

<KPj) = ~ , i = 1, 2, . . . , , </,(!/")= <MC/t n H , n=l,2,... . (47) 

Then 

^Ph^Pj^ ■ ■ ■ Pj n U^UP im ■ ■ ■ UP l2 UP n ) = 6 nm ^6 lin ■ ■ ■ 6 injn . (48) 
Proof. Put Qj = Pj - l/k, then 

Q t Qj = SijQi - k- l {Qi + Qj) + S^k- 1 - k~ 2 , 0(Q 3 -) = . (49) 
Hence the LHS of Q48JI is a linear combination of the terms of the form 

<l>(U* ni Q kl rf na --U mi Q kl U rn »--) , (50) 

which due to the freeness and Eqs. iHTlliSjl are all equal to zero except the terms which do not contain nontrivial 
powers of U and U'. This can happen, however, for n = m only. In this case we can easily prove H48|l by induction. 

The relation I|48J1 corresponds to the phenomenon observed in the numerical computations of St[X.,U] and 
Et [X, U, Vd\ and described by 14311 and 1441 . This justifies our hypothesis formulated above. In the case b) this 
hypothesis is not surprising due to the random choice of the partition and the generic relations between free random 
variables and random matrices. Similarly, we would expect the same phenomenon for the fixed partition and the 
random choice of the unitary matrix. On the other hand, for the case a) it seems to be a new characterization of 
chaotic quantum systems in terms of "quantum-probabilistic" relations between the dynamics and the measurement 
(coarse-graining) procedure. 



11 



8 Concluding remarks 



We have analyzed the decoherence in an open quantum system, the classical analogue of which is chaotic. The 
decoherence may be quantified by the rate of increase of the von Neumann entropy of the initially pure states. 
We have found an explicit upper bound for the rate of the von Neumann entropy given by the partial entropy 
used to define the ALF dynamical entropy. The later quantity is related to the Kolmogorov-Sinai entropy of the 
corresponding classical system. Hence our findings allow us to established a further relation between the speed of 
decoherence in open quantum systems and the degree of classical chaos. 

Such a relation, demonstrated in several earlier works |?fl ll()llTT| holds if some additional assumptions concerning 
the coupling of the system investigated with an environment (the measurement process) are made. In particular, 
we proposed to study the scheme of random measurements, in which the usual Kraus operators Xi, which represent 
projectors on some well defined fragments of classical phase space, are replaced by operators obtained by random 
matrices, X- — VXiV' . In such a case the rate of von Neumann entropy becomes maximal (with probability one, 
with respect to the choice of random matrix V). Thus the decoherence depends only on the kind of the measurement 
performed (the number of the Kraus operators or the dimensionality of the system), and is independent of the 
quantum unitary dynamics U, and of the degree of chaos (Lyapunov exponent, KS dynamical entropy) of the 
corresponding classical system. 

From a practical point of view it is therefore natural to ask, for which class of measurement procedures the 
relation between classical chaos and the degree of quantum decoherence is still valid. Although we are not in 
position to formulate mathematically rigorous sufficient conditions, which would imply such a relation, our numerical 
evidence allows us to advance the following conjecture. The interaction with an environment induces decoherence 
related to the degree of the classical chaos, if the measurement (Kraus) operators have a well defined classical limit. 
In other words, the coherent states representation of each of Kraus operators needs to be well localized in certain 
fragments of the classical phase space. 

More formally, the maximal entropy growth, and its independence of the unitary dynamics, may be analytically 
derived from an assumption that unitary operator U and an arbitrary combination of the projector operators Pi 
are free independent. Obviously this statement is of a statistical nature, and does not allow one to draw rigorous 
conclusion for a concrete set of projection and evolution operators. The free random variables approach concerns 
entire ensembles of operators and enables us to formulate exact statements concerning the decoherence rate in the 
limit of large Hilbert space dimension. Nevertheless, for practical purposes one may choose a set of arbitrary test 
states 4> and check whether property 1481 is approximately fulfilled for the analyzed unitary map U and measurement 
P. It is worth to emphasize that the free-independence condition can be formulated as a condition for a pair of 
genuinely quantum objects - an observable and an unitary quantum map - without any reference to the classical 
notions. Therefore the idea of quantum chaos may be extended to systems without obvious classical counterparts 
or to dynamics which do not satisfy standard assumptions concerning the spectral fluctuations. 
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acknowledged. 

References 

[1] A. Katok and B. Hasselblatt, Introduction to the modern theory of dynamical systems, Cambridge Univ. Press, 
Cambridge UK 1996. 

[2] M. Ohya and D. Petz, Quantum entropy and its use, Springer, Berlin 1993. 

[3] R. Alicki and M. Fannes, Quantum dynamical systems, Oxford Univ. Press, 2001. 

[4] W. Slomczyhski and K. Zyczkowski, J. Math. Phys. 35, 5674 (1994). 

[5] A. Connes, H. Narnhofer and W. Thirring, Commun. Math. Phys. 112, 691 (1987). 

[6] R. Alicki and M. Fannes, Lett. Math. Phys. 32, 75 (1994). 

[7] W. Slomczyhski and K. Zyczkowski, Phys. Rev. Lett. 80, 1880 (1999). 



[8] F. Benatti, V. Cappellini, M. De Cock, M. Fannes and D. Vanpeteghem, preprint quant-ph/0308069 



12 



[9] W. H. Zurek and J. P. Paz, Phys. Rev. Lett. 72, 2508 (1994). 
[10] D. Monteoliva and J. P. Paz, Phys. Rev. Lett. 85, 3373 (2000). 
[11] P. Bianucci, J. P. Paz, and M. Saraceno, Phys. Rev. E 65, 046226 (2002). 



[12] I. Garcia-Mata, M. Saraceno and M. E. Spina, preprint nlin. CD/0301025 

[13] F. Benatti, V. Cappellini and F. Zertuche, preprint |math-ph/ 0308033 

[14] R. Alicki, D. Makowiec and W. Miklaszewski, Phys. Rev. Lett. 77, 838 (1996). 

[15] N. L. Balazs and A. Voros, Ann. Phys. (N.Y.) 190, 1 (1989). 

[16] M. Saraceno, Ann. Phys. (N.Y.) 199, 37 (1990). 

[17] M. Saraceno and A. Voros, Chaos 2, 99 (1992); Physica D 79, 206 (1994). 

[18] S. De Bievre, M. Degli Esposti and R. Giachetti, Com. Math. Phys. 176, 73 (1996). 

[19] R. O. Vallejos and M. Saraceno, J. Phys. A: Math. Gen. 32, 7273 (1999). 

[20] A. Lozinski, P. Pakonski and K. Zyczkowski, Phys. Rev. E 66, 065201 (R) (2002). 

[21] R. F. Streater, Statistical dynamics, Imperial College Press, London 1995. 

[22] D. W. Robinson and O. Bratteli, Operator algebras and quantum statistical mechanics 2, Springer Verlag, 
Berlin 1997. 

[23] D. Voiculescu, Invent. Math. 104, 201 (1991). 

[24] D. Voiculescu, K. Dykema and A. Nice, Free random variables, CRM Monograph Series vol.1, American Math- 
ematical Society 1992. 



13 



